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TECHNICAL  REPORT  SUMMARY 

Objective 

The  purpose  of  this  study  is  to  mathematically  model  the  load-trunsfer 
characteristics  between  rock  and  selected  support  systems,  and  to  investigate 
the  response  to  varying  rock-behavioral  characteristics. 

General  Approach  and  Technical  Results 

Three  geotechnical  engineering  of  commercial  value  were  investigated. 

D  •<  jiled  reports  are  given  in  Appendices  A,  B,  and  C.  Roof  and  wall  movements 
a  an  underground  power  plant  were  calculated  and  composed  with  field  data. 
In-si tu  rock  mass  properties  were  determined  from  measured  roof  deflections, 
and  axial  and  berfding  moments  in  a  steel  support  were  calculated  to  determine 
the  applicability  of  two-dimensional  finite  element  computer  codes  to  that 
type  of  problem..  With  the  experience  gained  from  these  studies  a  computer 
code  was  developed  which  consolidated  all  the  existing  and  proven  finite 
element  technology  for  two-  and  three-dimensional  computer  cooes  into  one 
"user  oriented"  computer  code. 

Technical  Problems 

In  setting  up  and  executing  these  problems,  it  was  determined  that 
modeling  of  the  structural  features  of  the  rock  using  finite-element  computer 
codes  has  out-stripped  the  availability  of  rock  property  data. 

POD  Implications 
None . 

Implications  for  Further  Research 

Additional  research  should  be  made  in  the  field  to  determine  in-situ 
conditions  of  rock  mass  behavior  and  the  state  of  stress.  This  information 
can  then  be  used  in  existing  finite  clement  computer  codes. 

Special  Comments 


None. 
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THEORETICAL  INVESTIGATION  OF  ROCK  AND  SUPPORT  INTERACTION 

The  purpose  of  this  study  Is  to  mathematically  model  the  load-t' unsfer 
characteristics  between  rock  and  selected  support  systems,  and  to  Investigate 
the  response  of  these  structures  to  varying  rock-behavioral  characteristics. 

The  schedule  and  tasks  required  to  accomplish  this  mission  is  shown  in  figure  1. 
Tasks  undertaken  include:  (1)  determine  the  state-of-the-art  of  rock — support 
interaction  models;  (2)  determine  the  availability  of  engineering  parameters; 

(3)  develop  a  weighting  process  for  the  parameters;  (*0  obtain  and  modify  com¬ 
puter  codes;  and  (5)  perform  parameters  study.  Through  contract  monitoring 
(task  6),  computer  codes  were  developed  to  model  rock  mass  and  ground  support 
systems.  This  technology  was  then  applied  to  problems  of  commercial  value. 

Manpower  shortages,  employment  freeze,  and  the  reduction-in-force  pre¬ 
vented  full  accomplishment  of  the  tasks  delineated. 

Tasks  1  through  were  completed.  Wi tli  respect  to  tasks  1  through  A,  the 
following  conclusions  were  reached. 

Modeling  of  the  structural  features  of  rock  using  finite-element  computer 
codes  has  out-stripped  the  availability  of  rock  property  data.  Because  of 
the  complexity  of  rock  mass,  it  is  doubtful  that  detailed  finite-element  calcu¬ 
lations  will  ever  be  practicable  for  routine  design  of  underground  openings. 
Generalized  models  incorporating  only  the  salient  features  of  the  rock  mass 
could  be  used  for  phenominological  studies.  Finite-element  techniques  can 
be  helpful  in  computing  the  state  of  stress  for  structures  whose  behavior  is 
unsatisfactory  when  known  displacements  serve  as  a  check  on  the  values  of 
physical  properties  chosen.  Definite  need  exists  for  artificial  support  property 
data;  and  capability  of  3-D  models  must  be  improved  and  extended  to  increase 
the  scope  of  practical  problems  that  can  be  analyzed. 
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A  study  also  was  conducted  on  circular  and  square  openings  in  layered  rock. 

The  magnitude  of  the  field  stress  used  was  that  expected  in  deep  mining  in  the 
Coeur  d'Alene  mining  district  as  no  agreement  on  the  magnitude  and  direction  of  the 
stress  field  could  ba  found  in  literature.  The  physical  properties,  chosen  to  give 
an  in  situ  overall  elastic  modulus  of  2.5  x  10^,  and  the  allowable  compression  and 
tension  values  were  based  on  physical  property  tests  made  at  the  University  of  Idaho. 
Trend  of  th'  results  was  as  expected;  that  is  with  stiffer  material  there  is  less 
displacement;  there  is  a  preferred  orientation  of  the  opening  with  respect  to  the 
stress  field;  and  the  steeper  the  Mohr's  envelope,  the  less  plastic  deformation. 
Significance  of  this  study  is  that  under  similar  conditions  of  stress  and  rock 
properties,  rock  displacements  in  the  district  can  be  18  inches  while  the  elastic- 
plastic  calculations  give  only  3  inches — some  phenomena  other  than  plastic  flow  is 
the  primary  cause  of  the  rock  movement  in  field  situations.  Two  manuscripts  were 
written  describing  the  finite  element  analysis  of  rock  movement  in  an  underground 
powerplant  and  roof  movements  due  to  increasing  roof  span.  These  documents  are 
in  Appendix  A  and  Appendix  B.  Task  5  (to  begin  in  FY73)  was  incorporated  into  a 
statement  of  work  change  and  included  as  amodif ication  to  Ohio  State  University 
Contract  No.  H02i0017-^.  The  results  of  the  task  5  parametric  investigation 
showed  that: 

The  load  transfer  to  the  structural  support  system  is  dependent  upon  the 
tendency  of  the  rock  to  deform.  Higher  modulus  implies  less  deformation  for  the 
same  rock  load  and  consequently  less  load  transferred  to  the  steel  member. 

Higher  values,  tending  toward  .5,  of  Foisson's  ratio  are  associated  with 
redistribution  of  stress  in  the  rock.  This  stress  redistribution  is  reflected 
in  a  more  uniform  stress  distribution  in  the  steel  support. 

1/  Contract  H0210017  "Stresses,  Deformations,  and  Progressive  Failure  on  Non- 
Homogeneous  Rock"  with  Ohio  State  University  was  awarded  on  February  1,  1971. 

The  contract  was  later  amended  to  include  the  modification  mentioned  (Task  5). 
Information  on  obtaining  reports  and  magnetic  tape  containing  this  compuLer 
program  is  shown  on  page  55. 
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Sparce  blocking  increases  bending  moments  in  the  support  structure.  The 
total  load  transferred  to  the  support  structure  decreases  as  the  amount  of 
blocking  between  the  steel  set  ar.-'  rock  increases.  The  wood  blocking  acts 
much  like  a  "back  packing"  material. 

A  complete  manuscript  describing  the  details  of  this  work  is  given  in 

Appendix  C. 

As  a  consequence  of  this  work,  it  was  recommended  and  accepted  that  a  con¬ 
tract  be  funded  which  would  consolidate  all  the  existing  and  proven  finite  element 
technology  for  two-  and  three-dimensional  computer  codes  into  one  user  oriented 
computer  code.  This  contract  was  funded  (ARPA  contract  H0220035)  and  completed. 
Appendix  D  is  a  descriptive  summary  of  the  computer  code  capability.  Subsequently, 
an  era  of  computer  code  application  is  anticipated.  From  application,  we  hope  to 
discover  the  deficiencies  within  the  state-of-the-art  and  identify  where  fundamental 
research  is  needed  and  make  appropriate  recommendations  for  future  research. 
Manuscripts: 

"Finite  Element  Analysis  of  Excavations  in  Rock,"  by  Chin-Yung  Chang, 

Keshavan  Nair,  and  William  J.  Karwoski.  Oral  presentation  and  subsequent 
publication  in  the  proceedings  of  Corps  of  Engineers,  Waterways  Experiment 
Station  Symposium  on  "Applications  of  the  Finite-Element  Method  to  Problems 
in  Geotechnical  Engineering",  held  May  l-*t  in  Vicksburg,  Mississippi. 

"Roof  Movement  Due  to  Increasing  Roof  Span  in  an  Underground  Opening, 
by  William  J.  Karwoski,  Ranbir  S.  Sandhu,  R.  Singh. 

Oral  presentation  and  subsequent  publication  of  abstract  in  the  proceed¬ 
ings  of  the  l*ith.  Annual  Rock  Mechanics  Symposium,  held  June  12-1*1,  in  Univer¬ 
sity  Park,  Pennsylvania. 


"A  Parametric  Study  of  Stresses  in  Steel  Support  for  a  Tunnel,"  by  R^nbir 
S.  Sandhu,  and  William  J.  Karwoski .  Oral  presentation  and  subsequent  publication 
in  proceedings  of  the  15th  Annual  Rock  Mechanics  Symposium,  held  September  17- 
1?  at  the  South  Dakota  School  of  Mines  and  Technology,  Rapid  City,  South  Dakota. 
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The.  manuscript  of  "Finite  Element  Analysis  of  Excavating  in  Rock" 
presented  by  the  Symposium  on  Applications  of  The  Finite  Element  Method 
in  Geotechnical  Engineering  held  at  the  U.S.  Army  Corps  of  Engineers, 
Waterways  Experiment  Station,  Vicksburg,  MS,  May  1-4,  1972. 

This  presentation  included  results  of  research  performed  under 
contract  H02 10046  "A  Theoretical  Method  for  Evaluating  Stability  of 
Openings  in  Rock"  with  Woodward-Lundgren  &  Assoc.,  Oakland,  Calif., 

Dr.  Keshavan  Nair,  Principal  Investigator.  Technical  reports  from  this 
contract  and  the  extension  of  this  work  (contract  110220038  with  the 
same  organization)  are  available  from  the  National  Technical  Information 
Service,  U.S.  Department  of  Commerce,  5285  Port  Royal  Road,  Springfield, 
Va.  22151.  Following  is  the  information  necessary  for  ordering  these 
reports . 
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Final  Technical  Report 
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FINITE  ELEMENT  ANALYSIS  OF  EXCAVATIONS  IN  ROCK 

By 

Chin-Yung  Chang^,  Keshavan  Nair^  and  William  J.  Karwoski^ 

SUMMARY 

A  general  finite  element  computer  program  for  the  analysis  of 
plane  problems  which  .incorporates  no  tension,  joint  perturbation 
and  elasto-plastic  behavior  was  developed.  The  -esults  of  anal¬ 
yses  were  compared  with  actual  performance  for  the  Morrow  Point 
Powerplant  excavation  and  a  laboratory  model  study  of  an  opening 
in  a  rock-like  material. 

The  comparison  indicated  that  the  developed  analytical  techniques 
had  the  capability  of  providing  a  good  qualitative  prediction  of 
performance.  Quantitative  prediction  was  also  considered  to  be 
acceptable  in  view  of  the  input  information  and  idealization. 


1.  Research  Engineer,  Woodward-Lundgren  5  Associates,  Oakland, 
California . 

2.  Vice-President,  Woodward-Lundgren  8  Associates,  Oakland, 
California 

3.  Mining  Engineer,  Bureau  of  Mines,  Spokane,  Washington. 
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mass  due  to  yielding  have  beer  presented  by  Reyes  and  Deere 
(1966).  Preliminary  analysis  has  indicated  that  these  new 
methods  of  analysis  using  finite  element  techniques  have  the 
potential  of  predicting  performance  with  improved  accuracy. 

There  has  been  very  limited  verification  cf  these  techniques 
on  the  basis  of  comparison  v.’ith  measured  field  performance. 
Without  field  verification,  the  use  of  these  new  techniques 
in  the  design  of  excavations  in  rock  will  remain  limited. 
Therefore,  to  develop  a  theoretically  sound  method  for  de¬ 
signing  excavations  in  rock,  which  will  be  used  in  practice, 
an  essential  step  is  to  verify  the  ability  of  the  available 
methods  of  stress  analysis  in  predicting  the  behavior  of  rock 
masses . 

The  study  reported  herein  can  be  divided  into  two  phases: 

(1)  Development  of  a  general  analytical  procedure  (i.e.,  a 
finite  element  computer  program  with  wide  capabilities)  for 
determining  the  mechanical  state  in  a  rock  mass,  and  (2)  Anal¬ 
ysis  of  case  histories  to  compare  predicted  and  observed  values. 

DEVELOPMENT  OF  A  GENERAL  COMPUTER  PROGRAM 

In  order  to  develop  a  single  general  computer  program,  it  was 
necessary  to  develop  a  consistent  computational  technique  to 
model  the  different  aspects  of  rock  behavior.  On  the  basis  of 
the  review  of  the  available  techniques,  it  was  concluded  that 
the  "initial  stress"  (stress  transfer)  technique  presented  by 
Zienkiewicz  and  his  co-workers  would  provide  a  consistent 
approach  in  the  development  of  a  finite  element  computer  program 
for  t ime - independent  plane  problems  which  included  the  following 
capabilities:  (1)  No  Tension  Analysis,  (2)  Joint  Perturbation 

Analysis,  and  (3)  Elasto-Plastic  Analysis.  The  major  reason  for 
selecting  this  approach  was  the  computational  advantage  that  re¬ 
sults  from  using  the  initial  elastic  stiffness  at  every  stage  of 
the  solution  process.  The  essential  concepts  used  to  include  the 
above  listed  rock  characteristics  are  discussed  in  the  subsequent 
paragraphs. 
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I .  No  Tension  Analysis 

The  basic  concept  used  in  the  combined  computer  program  for 
the  no  tension  analysis  is  similar  to  that  developed  by  Zien- 
kiewicz,  et  al  (1968).  The  major  steps  in  the  analysis  used 
can  be  summarized  as  follows: 

1.  Assign  initial  stresses  to  the  rock  mass,  and  cal¬ 
culate  the  boundary  loads  required  on  the  cavity 
face  to  simulate  the  creation  of  the  opening  and 
other  structural  loads  applied  to  the  system. 

2.  Analyze  the  problem  as  a  linear  elastic  problem. 

Add  the  induced  changes  in  stress  to  the  initial 
stresses  and  compute  the  principal  stresses. 

3.  Determine  those  elements  in  which  tension  exists. 

If  the  material  is  assumed  incapable  of  sustaining 
any  tension,  or  if  the  tensile  stres~  exceeds  the 
tensile  strength,  the  excess  tensile  stresses  have 
to  be  eliminated.  This  is  done  by  applying  nodal 
point  forces  calculated  to  eliminate  tensile  stresses. 


4.  The  elastic  analysis  is  repeated  for  the  calculated 
equivalent  nodal  point  forces;  stresses  are  determined 
in  the  elements  The  check  for  tensile  stresses  is 
repeated . 

5.  If,  at  the  end  of  step  (4),  principal  tensile  stresses 
are  still  present,  steps  (3)  and  (4)  are  repeated  until 
there  is  no  appreciable  difference  in  magnitude  and 
distribution  of  stresses  with  further  iterations. 


In  Step  (3)  when  the  linear  elastic  solution  indicates  that  the 
rock  is  subjected  to  tensile  stresses  greater  than  the  tensile 
strength,  the  rock  is  assumed  to  be  fractured  and  incapable  of 
transferring  stresses  between  two  orthogonal  directions.  To 


include  this  effect,  a  correction  has  to  be  made  on  the  stress 

before  the  stress  transfer  process  is  performed.  In  the  case 

where  both  principal  ten? ile  stresses  o'  and  a  '  are  to  be 

transferred,  the  'corrected'  stresses  are  given  by  (Chang  anc'  Nair 
1972) 

°x  *=  (1  +  v2  +  .  .)  ox'  +  v  (3  +  v2  +  .  .  .)  o  ' 

(1) 

°y  c  (1  +  v2  +  ...)  oy'  +  v  (1  +  v2  +  ...)  ox' 

II  •  Joint  Perturbation  Analysis 

The  stiffness  matrix  of  one-dimeni> Awnal  joints  is  formulated 
according  to  the  procedure  developed  by  Goodman,  Taylor,  and 
Brekke  (1968).  The  one-dimensional  joint  with  its  local  coor¬ 
dinate  system  and  the  sign  convention  of  its  normal  and  shear 
stresses  is  illustrated  in  Figure  1. 

The  shear  strength  of  a  joint  ma/  be  expressed  by 

Tf  «=  c  +  oN  tan  <f,e  (2) 

where 

t£  -  shear  strength  of  the  joint 
c  *=  'cohesion  along  the  joint 
°N  “  normal  stress  across  the  joint 
4*e  c  effective  friction  angle  of  the  joint  surface 
(Patton  1966) 

If  tne  normal  stress  across  the  joint  is  tensile,  it  is  assumed 
that  the  joint  is  incapable  of  resisting  any  shear  stress,  i.e.  , 
it  has  no  strength.  A  flow  diagram  illustrating  the  stress  trans¬ 
fer  technique  for  the  joint  perturbation  analysis  is  shown  in 
Figure  2. 


Ill*  Elasto-Plastic  Analysis 

In  developing  an  eiusto-piastic  analysis,  it  is  necessary  to 
define  a  yield  function  and  the  stress-strain  relations  before 


and  after  yield.  Prior  to  yield,  it  is  assumed  that  linear 
elastic  stress-strain  relations  are  applicable. 


Yield  Fun;tion 

The  yield  function  utilized  is  a  generalization  of  the  Mohr- 
Coulomb  hypothesis  suggested  by  Drucker  and  Prager  (1952). 

The  yield  function  is  represented  by  the  following  equation: 

f  **  a J i  +  JTT  =  k  (3) 


where : 

f  and  k  =  material  constants 

Ji  *=  first  stress  invarient 

J2  *  second  invarient  of  stress  deviation. 

pt*  =  material  constant 

The  yield  surface  expressed  by  Equation  (3)  for  a>0  is  a  right 
circular  cone  with  its  axis  equally  inclined  to  the  coordinate 
axis . 


In  the  case  of  plane  strain,  Drucker  and  Prager  (1952)  have 
shown  that 

a  -  - -  (4) 

(9  +  12  tan2  4 >)'s 

and 

k  -  - — - -  (5) 

(9  +  12  tan2  ({O'* 

where : 

c  ■  the  cohesion  of  the  material 

<(>  ■  angle  of  internal  friction  of  the  material. 

Stress-Strain  Relations 

During  an  infinitesimal  increment  of  stress,  changes  in  strain 
are  assumed  to  be  composed  of  clastic  and  plastic  parts  if  the 
element  is  in  yield,  i.e. 
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(Ac)  e  (Ac)  +  {Ac}  (6) 

e  P 

The  elastic  strain  increments  are  related  to  stress  increments 
by  the  generalized  Hooke's  law  as 

{Aa}  =  [D]  {Ae)e  (7) 

[D]is  the  linear  isotropic  elastic  coefficient  matrix. 

Utilizing  the  Drucker  Prager  criteria  (equation  3)  ,  Reyes  (1966) 
developed  elasto-plastic  stress-strain  relations.  These  rela¬ 
tions  can  be  expressed  as  follows: 

(4o)  -  [D]e  p  (4e)p  '  (8) 


(1  +  9ct: 


[•-*]  [“ 


3Ka  J] 

*TT  ‘  77 


(1  +  9a2  £) 


3\>Kk _ 

V5  (1  *  9c!  £) 


(10) 


l-  =  - 

2 J 2  (1  +  9a2  £) 


G  (shear  modulus)  = 


K  (bulk  modulus)  = 


(11) 


Since  the  elas to-plastic  stress  -  strain  relation  is  a  function 

of  the  stress  state,  it  is  necessary  to  keep  a  record  of  the 

axial  stress  a  .  In  the  elastic  range  under  the  plane  strain 
z 

conditions,  the  following  relationship  exists  between  the  change 
in  axial  stress  (Aa  )  and  the  changes  in  stress  in  the  x  -  y  plane. 


Aa  =  v(Aox  +  Aa2) 

id 


(12) 


where  Ao x  and  A o2  are  changes  in  stress  in  two  principal  stress 
directions  in  x  -  y  plane.  In  the  plastic  range,  Drucker  and 
Prager  (1952)  have  shown  that  the  change  in  axial  stress  is  given 
by  : 


A°z  =  J  (Aoi  +  Aa2)  -  j  (Aai  -  Ao2)  Sin  $ 


(13) 


To  conduct  an  incremental  elasto-plastic  analysis,  the  load  is 
applied  in  a  series  of  increments.  A  flow  diagram  illustrating 
the  computational  procedure  using  the  stress  transfer  technia.ue 
is  shown  in  Figure  3. 
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FIG.  3  -'STRESS  TRANSFER'  TECHNIQUE  FOR  ELASTO-PLASi If  ANALYSIS 
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A  macro  flow  chart  indicating  the  organization  of  the  complete 
program  is  presented  in  Figure  4. 

ILLUSTRATIVE  EXAMPLES 

The  following  three  problems  were  solved  to  verify  the  program 
an  illustrate  its  use: 

I.  Elasto-plast ic  analysis  of  a  thick-walled  circular  tube 
with  the  Von  Mises  yield  criterion.  A  closed  form  solu¬ 
tion  is  available  for  this  case  for  verification,  Prager 
and  Hodge  (1951 ) . 

II.  Elasto-plastic  analysis  of  a  circular  opening  with  the 
generalized  Mohr-Coulomb  yield  criterion.  The  results 
are  compared  with  those  obtained  by  Reyes  (1966)  and 
Baker,  et  al  (1969) . 

III.  Combined  no  tension,  joint  perturbation  and  elasto- 
plastic  analysis  of  a  rectangular  underground  opening 
to  demonstrate  the  usage  of  the  combined  computer 
program. 

The  results  of  these  analyses  have  been  compared  with  closed 
form  solutions  and  the  solutions  obtained  by  other  researchers. 
Typical  results  are  presented  in  the  Appendix. 

ANALYSIS  OF  LABORATORY  MODEL 
Description  of  Model  Study 

Heuer  and  Hendron  (1971)  conducted  a  series  of  model  tests  under 
plane  strain  conditions  on  unlined  openings  in  a  rock-like  mate¬ 
rial.  The  geometry  of  the  model  block  which  consisted  of  two 
halves  assembled  together  is  shown  in  Figure  5a.  The  behavior  of 
the  model  was  monitored  by  internal  foil  strain  gage  rosettes 
installed  on  the  block  interface  at 'various  radii.  Circumferential 
strain  gages  were  also  installed  on  the  wall  of  the  opening.  Two 
sets  of  diametrical  extensometers  were  also  used  to  measure  the 


MAIN  PROGRAM 


SIMPLIFIED  FLOW  DIAGRAM  SHOWING  SEQUENCE  OF  OPERATION  OF  ALL  SUBROUTINES 
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closure  of  the  opening.  The  loads  applied  to  the  model  con¬ 
sisted  of  horizontal  and  vertical  stresses  applied  at  the 
boundary  as  shown  in  Figure  5a.  Tests  were  conducted  under 
three  ratios  of  horizontal  to  vertical  stress,  (a^/^v) ,  1.0, 

.25  and  .75.  Only  the  tests  with  oh/oy  -  .25  are  analyzed  in 
this  stud/. 

Based  on  the  results  obtained  from  triaxial  tests  on  the  mate¬ 
rial  utilized  for  constructing  the  model  ( Heuer  and  Hendron  , 

1971)  the  material  constants,  (a,  k)  used  in  Equation  (3)  for 
the  yield  function  were  calculated  to  be  0.262  and  189  psi, 
respectively.  The  tensile  strength  of  the  material  was  deter¬ 
mined  to  be  35  psi  as  compared  with  the  unconfined  compressive 
strength  of  600  psi.  The  axial  s tress -strain  curves  and  the 
value  of  Poisson's  ratio  were  found  to  be  functions  of  the 
principal  stress  difference.  The  stress-strain  curves  indicate 
that  a  linear  approximation  is  valid  up  to  a  principal  strep* 
difference  (oj  -  03)  of  600  psi.  A  linear  isotropic  elastic 
characterization  was  utilized  for  the  analysis. 

During  testing  at  a^/ov  =  1/4 ,  fractured  wedges  occurred  at 
the  snringlines  when  the  applied  load  was  incremented  from 
°v  =  690  Psi  to  °v  =  ^65  psi,  Figure  5b.  As  the  stress  level 
was  increased,  more  fractures  formed  and  the  wedges  moved  into 
the  opening.  Another  significant  set  of  fractures  were  formed 
extending  away  from  the  opening  back  into  the  model.  These 
fractures  appear  to  be  extensions  cf  the  fractures  which  formed 
the  wedges  at  the  springline,  Figure  5b.  At  the  crown  and  invert 
high  circumferential  tensile  strains  were  measured  at  applied 
stress  levels  above  oy  -  400  psi.  Heuer  and  Hendron  (1971) 
hypothesized  that  tension  cracks  might  have  formed  during  testing 
However,  no  cracks  were  visible  after  testing.  It  was  assumed 
that  the  cracks  were  open  during  the  test  and  then  closed  upon 
removal  of  load. 


Analysis  of  Model  Study 

Idealization  of  the  Model  -  Because  of  the  symmetry  of  the 
model  with  respect  to  its  axes,  as  shown  in  Figure  5a,  it  was 
only  necessary  to  analyze  a  quadrant  of  the  cross-section. 

The  finite  element  idealization  and  the  elastic  parameters 
used  in  the  analysis  are  shown  in  Figure  6. 


The  incremental  analysis  was  performed  to  simulate  the  actual 
loading  path.  The  incremental  loads  applied  are  shown  in 
Figure  6.  Two  cases  were  analyzed  under  the  assumption  that 
the  material  behaved  as  an  elasto-plastic  material  and  had  a 
limited  tensile  strength  of  35  psi.  The  first  case  (Case  1) 
assumed  that  the  block  acted  as  a  continuum.  The  second  case 
(Case  2)  assumed  that  a  fractured  wedge  formed  at  the  spring¬ 
line  during  testing  and  moved  into  the  opening. 


Presentation  and  Pircussicr.  of  P.ccults 


Figure  7  illustrates  the  development  of  plastic  and  tensile 
regions  at  three  levels  of  the  axial  stress  400,  600,  and 
800  psi.  The  observed  fractures  are  also  shown  in  Figure  7. 

It  may  be  noted  in  Figure  7  that  the  development  of  the  plastic 
region  calculated  at  oy  *=  600  psi  is  similar  to  the  shear- 
fractured  wedges  developed  at  the  springline.  At  ov  =  800  psi, 
a  large  plastic  region  developed  back  into  the  model.  The  loca¬ 
tion  of  the  actual  shear  fracture  is  shown  to  be  along  the  edge 
of  the  plastic  region.  It  was  observed  in  the  model  test  that 
a  region  of  tensile  stress  developed  at  the  crown  and  invert. 
However,  there  was  little  failure  in  this  zone  indicating  that 
the  tensile  stresses  were  below  the  tensile  strength  of  the 
material.  An  elasto-plastic  analysis  indicates  that  a  small 
zone  of  tension  cracks  would  develop  if  the  material  exhibits 
elasto-plastic  behavior.  The  results  of  the  elasto-plastic 
analysis  appear  to  agree  with  observed  behavior  in  the  model  test. 
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In  order  to  simulate  the  propagation  of  fractures  and  their 
effects  on  the  behavior  of  the  model,  an  analysis  was  performed 
(Case  2)  after  removing  the  wedge  that  formed  under  a  vertical 
stress  of  600  psi.  Under  these  conditions,  the  plastic  zone 
extended  behind  the  wedge  and  propagated  into  the  model  block 
as  shown  in  Figure  8.  The  direction  of  the  propagation  of  the 
plastic  zone  and  the  actual  propagation  of  the  fracture  as  the 
loads  increased  were  similar.  The  extend  of  the  plastic  zone 
as  the  vertical  stress  was  increased  to  800  psi  is  also  shown 
in  Figure  8.  Also  shown  is  the  fracture  which  developed  during 
the  test  as  illustrated  in  Figure  5b.  It  may  be  noted  that  the 
shear  fracture  observed  lies  in  the  middle  or  the  plastic  zone 
indicating  the  analysis  can  provide  a  good  indication  with  respect 
to  the  development  of  the  critical  zone  in  the  vicinity  of  the 
opening . 

If  it  is  recognized  that  the  yield  function  can  also  represent 
a  limiting  stress  state  for  fracture  to  occur,  then  the  results 
of  the  elasto-plastic  analysis  are  in  good  agreement  with  observed 
performance . 

In  order  to  compare  with  the  strains  measured  in  the  model  tests, 
element  strains  were  computed  from  the  nodal  point  disr lacements . 
These  strains  were  assumed  to  occur  at  the  centroid  of  the  element. 
Computed  vertical  and  horizontal  strains  along  the  line  (0  -  86.25°) 
close  to  the  horizontal  axis  (6  =  90°)  for  Cases  1  and  2  are  shown 
in  Figures  9  and  10,  respectively.  For  the  purpose  of  comparison, 
the  measured  tangential  and  radial  strains  along  the  horizontal 
axis  are  also  plotted  in  Figures  9  and  10. 

' 

Because  of  the  lack  of  reliability  of  the  test  data  obtained  from 
the  electrical  strain  gages  as  evidenced  by  the  relatively  large 
scatter  in  the  measurements,  Heuer  and  Hendron  (1971)  concluded 
that  the  strain  measurements  would  only  give  a  qualitative  indi¬ 
cation  of  the  response  of  the  model.  Comparison  of  the  measured 
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DISTRIBUTION  OF  HORIZONTAL  STRAINS  ALONG  6  =  86.25 
ANALYSIS  OF  MODEL  TESTS 


and  calculated  behavior  shows  that ,  at  the  low  stress  levels, 
they  are  in  good  agreement.  However,  as  the  stress  level 
increased  and  the  fractured  wedges  formed,  the  strains  in¬ 
creased  to  a  level  much  greater  than  that  predicted  from  the 
Case  1  analysis  in  which  the  model  block  was  assumed  to  be 
a  continuum.  The  results  of  the  Case  2  analysis  indicated  a 
strain  pattern  and  magnitude  in  approximate  agreement  with 
measured  values.  This  indicates  that  the  formation  of  fractured 
wedges  at  the  springline  has  a  great  effect  on  the  behavior  of 
the  model  block.  Discrepancies  between  the  measured  strains 
and  the  strains  computed  by  Case  2  analysis  at  the  higher 
stress  levels  may,  in  large  part,  be  attributed  to  the  for¬ 
mation  of  the  second  set  of  fractures  extending  into  the 
model  block  away  from  the  opening. 

Analysis  of  the  model  tests  indicates  -that  the  elasto-plastic 
analysis  may  provide  valuable  information  for  evaluating  the 
critical  zone  in  an  intact  rock  surrounding  an  opening  Although 
the  analysis  cannot  simulate  accurately  the  behavior  of  the  model 
after  the  formation  of  shear  fractures  surrounding  the  opening, 
some  indication  can  be  obtained  with  regard  to  propagation  of 
fractured  zones. 


The  Morrow  Point  Underground  Powerplant  was  constructed  by  the 
U.S.  Bureau  of  Reclamation  on  the  Gunnison  River  some  20  miles 
east  of  Montrose,  Colorado  (Dodd,  1967  and  Brown,  et  al.  1971). 

The  powerplant  chamber  is  206  ft  long  and  57  ft  wide  with  a 
height  ranging  from  65  to  134  ft  and  is  located  about  400  ft 
below  the  ground  surface.  The  plan  of  the  powerplant  and  other 
adjacent  structures  is  shown  in  Figure  11.  A  cross-section  of 
the  powerplant  chamber  along  A-A1 (line  4+12  ft)  is  shown  in 
Figure  12.  It  may  be  noted  from  Figures  11  and  12  that  the  power- 
plant  chamber  is  situated  behind  a  steep  valley  rock  wall.  The 
crown  of  the  chamber  on  the  river  side  (b  line  wall,  Figure  12) 
is  about  200  ft  behind  the  ground  surface.  Details  with  regard 
to  the  geology  of  ;he  rocks  are  given  by  the  Bureau  of  Reclama¬ 
tion  (1965) . 


There  are  two  distinct  shear  zones  intersecting  the  powerplant 
area.  The  lower  zone  (shear  zone  A)  strikes  N  40°  W,  dips 
J2°E,  and  consists  of  a  one  to  five  ft  zone  of  fault  gouge 
and  fractured  rock.  The  upper  zone  (shear  zone  B)  has  a 
strike  of  N  ?0°  E,  a  dip  of  23°  E  and  consists  of  a  one  to 
three  ft  zone  of  fault  gouge  and  fractured  roc...  The  orien¬ 
tations  of  the  three  major  joint  sets  that  intersect  the 
chambers  are:  strike  N  63°  W  and  dip  82°  SW;  strike  N  36°  E, 

and  dip  80°  W;  strike  N-S  and  dip  43°  E. 

Significant  Features  of  the  Behavior  of  the  Powerplant  Excavation 
The  crown  of  the  chamber  was  first  excavated  during  a  three-month 
period  from  May  of  1965  through  July  of  1965.  The  excavation  of 
the  remaining  chamber  continued  from  August  of  1965  through  March 
of  1966.  When  the  excavation  reached  the  El.  6793  bench,  48  ft 

below  the  crown,  some  initial  movement  of  the  a-line  rock  wall 

occurred.  When  the  powerplant  was  excavated  to  El.  6748,  93  ft 

Kn  1  ru.»  4-  V>  n  J  J  i.  ->/*  1  r  - 1  A  - 
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was  measured  at  line  4+12  ft.  As  the  excavation  proceeded, 
the  a-'line  rock  wall  continuously  moved  inward  but  at  a  diminishing 
rate.  Significant  rock  movement  ceased  by  the  end  of  March,  1966 
(Brown,  Morgan  and  Dodd,  1971).  The  maximum  observed  inward 
movement  of  the  a-line  wall  at  line  4  +  12  ft  was  2.5  inches. 

Figure  13  shows  the  history  of  the  observed  rock  movements  at 
line  4  +  12  ft  reported  by  Brown,  et  al  (1971). 


Brown,  et  al.  (1971)  hypothesized  that  the  anomalous  rock  be¬ 
havior  along  the  a-line  rock  wall  might  involve  a  mass  movement 
of  a  rock  wedge  shown  in  Figure  14.  It  is  hypothesized  that 
shear  zones  A  and  B  intersect  at  an  average  distance  of  130  ft 
behind  the  a-line  wall,  forming  a  rock  wedge.  In  addition, 
two  planes  of  failure  must  exist  for  the  wedge  to  move  into 
the  opening.  One  is  the  upper  shear  failure  plane  with  an 
apparent  dip  of  17°  intersecting  the  chamber  near  the  spring¬ 
line  of  the  rock  arch.  The  other  is  essentially  vertical  and 


INDICATES  MOVEMENT  TOWARDS  CENTERLINE  OF  CHAMBER 
INDICATES  MOVEMENT  AWAY  FROM  CENTERLINE  OF  CHAMBER 
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DIMENSION  VIEW  OF  ROCK  WEDGE,  MORROW  POINT  POWERPLANT 
.  BROWN,  ET  AL.  1971) 
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intersects  the  chamber  between  line  2  +  15.5  ft  and  the  east 
end  wall  of  the  chamber.  A  cross-section  through  the  chamber 
at  line  4  +  12  ft  showing  the  rock  wedge  is  presented  in 
Figure  15. 


Idealization  of  the  Powerplant  Excavation 

a.  Finite  Element  Idealization  -  As  discussed  previously  the 
observed  maximum  rock  movements  occurred  at  line  4  +  12  ft. 
Therefore  a  cross-section  through  the  opening  at  this  location 
(Figure  12)  was  chosen  for  analysis.  The  essential  features 
in  the  idealization  are  the  steep  sloping  valley  ground  surface, 
shear  zones  A  and  B  and  two  other  incipient  failure  planes. 

The  steep  sloping  ground  surface  makes  the  estimation  of  the 
initial  in-situ  stresses  difficult.  As  described  previously, 
both  the  rock  wedge  and  the  movements  associated  with  it  are 
three-dimensional  in  nature.  However,  for  the  purposes  of 
analysis,  it  was  necessary  .o  idealize  it  as  a  two-dimensional 


plane  strain  problem.  Shear  zoned  A,  £  etiiu  the  other  twu  in¬ 
cipient  failure  planes  are  shown  as  distinct  discontinuities 
extending  in  the  direction  perpendicular  to  the  cross-section. 
These  discontinuities  are  idealized  as  Goodman's  one-dimensional 


joints.  Its  deformability  is  approximated  by  normal  and  shear 
stiffness . 


b.  Material  Properties  -  Extensive  testing  programs  were  con¬ 
ducted  by  the  Bureau  of  Reclamation  (1965)  to  determine  rock 
properties  at  the  Morrow  Point  Dam  Site.  These  included  direct 

shearing  and  sliding  friction  both  in  the  field  and  laboratory, 

*  * 

triaxial  testings  and  field  jacking  tests  of  foundation  rock. 

The  strength  and  elastic  properties  of  the  rock  above  and  below 
shear  zone  A  used  in  the  analysis  are  summarized  in  Table  1. 

As  described  previously,  both  shear  zones  A  and  B  were  idealized 
by  one-dimensional  joint  elements.  The  properties  of  these 
shear  zones  are  approximated  by  the  normal  and  shear  joint 
stiffnesses  (Goodman,  1969)  which  are  functions  of  normal  and 
tangential  deformability  of  the  shear  zones  and  their  geometries. 


FIG.  15  -  SCHEMATIC  DIAGRAM  OF  A-LINE  ROCK  WEDGE  MOVEMENT 
(AFTER  BROWN ,  ET  AL.  1971) 
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No  data  were  available  for  the  normal  and  shear  stiffnesses  of 
the  shear  zones  with  thickness  varying  from  1  to  5  ft.  There¬ 
fore,  two  sets  of  joint  stiffness  were  utilized  to  study  the 
influence  of  the  def ormability  of  the  shear  zones  on  the  be¬ 
havior  of  the  powerplant  excavation.  The  strength  parameters 
assumed  for  the  shear  zones  are  presented  in  Table  1. 

Because  of  the  limited  data  with  regard  to  the  magnitude  of 
the  in-situ  stresses  and  the  difficulties  associated  with 
estimating  their  values  in  the  rock  mass,  the  initial  state 
of  stress  used  in  the  analysis  was  obtained  by  a  gravity 
turn-on  analysis  of  the  rock  mass  without  the  opening. 

These  analyses  were  conducted  to  account  for  the  effects 
of  the  steep  valley  wall  located  in  the  vicinity  of  the 
powerplant  excavation.  For  these  analyses,  the  higher 
values  of  the  Poisson's  ratio  (=0.49)  was  assumed  for  the 
rock  to  simulate  the  likelihood  of  a  high  horizontal  stress. 

The  values  of  the  initial  stresses  which  have  to  be  applied 
at  excavated  boundary  to  simulate  the  excavation  were  obtained 
by  a  method  similar  to  that  suggested  by  Clough  and  Duncan 
(1969).  The  nodal  point  stresses  on  the  excavated  boundary 
were  estimated  from  the  stresses  in  the  surrounding  elements. 

Analysis  Procedures  . 

As  indicated  previously  two  cases  using  high  and  low  normal 
and  shear  joint  stiffnesses  for  the  shear  zones  were  analyzed. 
For  each  case,  the  initial  stresses  for  each  element  were  first 
obtained  by  performing  a  gravity  turn-on-analysis.  The  nodal 
stresses  were  then  determined  along  the  excavated  surface 
and  applied  at  the  excavated  boundary  to  simulate  the  excava¬ 
tion.  Very  small  values  of  the  elastic  constants  were  assigned 
to  those  elements  in  the  opening  to  simulate  the  cavity. 
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Presentation  and  Discussion  of  Results 

The  results  of  the  analyses  are  presented. in  Figures  16  and  17. 
Figure  16  illustrates  the  horizontal  displacements  fpr  points 
along  the  face  of  the  powerplant  chamber,  and  Figure  17  shows 
the  vertical  displacements.  The  observed  movements  on  a-  and 
b-line  walls  are  also  shown  in  Figures  16  and  17.  For  the  case 
with  high  joint  stiffnesses,  the  horizontal  inward  movement  at 
El.  6793  was  calculated  to  be  0.44  in.  on  the  a-line  wall  and 
0.41  in.  on  the  b-line  wall.  When  the  low  joint  stiffnesses 
were  used,  the  calculated  inward  movements  increased  to  1.35  in. 
on  the  a-line  wall  and  0.62  in.  on  the  b-line  wall. 

The  results  indicate  that  the  analysis  provided  a  good  quali¬ 
tative  estimate  of  rock  movements  at  the  Morrow  Point  Power- 
plant  excavation.  The  computed  deformations  are  of  the  same 
order  of  magnitude  as  those  observed.  The  difference  between 
the  computed  and  observed  deformation  is  significant  in  terms 
of  percentages  but  is  considered  satisfactory  in  terms  of 
numerical  values  when  it  is  recognized  that  the  following 
approximations  were  made  in  the  analysis  (1)  The  rock  wedge 
and  the  movements  associated  with  it  are  three-dimensional 
in  nature,  Figure  14.  A  two-dimensional  plane  strain  analysis 
would  tend  to  underestimate  the  movements.  (2)  The  initial 
state  of  stress  was  estimated  in  the  analyses.  It  is  quite 
possible  that  a  higher  initial  state  of  stress,  especially 
higher  horizontal  stresses  might  exist  in  the  rock  mass  which 
would  cause  higher  induced  stresses  and  movements.  (3)  No 
data  were  available  on  the  properties  of  the  shear  zone 
present  in  the  rock  mass.  It  was  necessary  to  assume  their 
values  in  the  analysis. 

It  should  also  be  recognized  that  the  powerplant  was  excavated 
in  stages  whereas  the  analysis  assumed  that  the  excavation  was 
created  instantaneously.  For  an  ideal  elastic  material  the 
assumption  would  not  result  in  any  error.  However,  in  those 
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cases  where  the  loading  path  has  an  influence  on  the  results 
the  effect  of  this  assumption  can  be  significant.  High  stress 
concentrations  during  the  excavation  process  could  cause  dis¬ 
placements  which  would  be  larger  than  those  computed  on  the 
assumption  of  an  instantaneous  creation  of  the  opening. 

CONCLUSIONS 

Computational  techniques  which  have  the  capability  of  including 
realistic  idealizations  of  in  situ  rock  conditions  can  provide 
a  good  qualitative  prediction  of  the  behavior  of  a  rock  mass 
in  the  vicinity  of  an  excavation.  Considering  the  availability 
and  accuracy  of  input  information  the  ability  of  the  developed 
computational  techniques  to  predict  quantitative  behavior  is 
considered  acceptable. 

It  was  found  that  the  effort  necessary -to  idealize  a  practical 
problem  for  analysis  is  considerable.  Since  real  problems,  in 
general,  include  certain  factors  which  cannot  be  modelled  accurately, 
it  may  be  necessary  to  utilize  more  than  one  idealization  to  obtain 
bounds  on  the  actual  behavior. 

Theoretical  concepts  and  computational  techniques  need  to  be 
developed  for  including  (i)  the  development  and  propagation  of 
fractures  in  a  rock  mass,  and  (ii)  the  influence  of  construction 
sequence  and  support  systems. 

It  is  recommended  that  a  major  effort  should  be  directed  towards 
the  study  of  case  histories  for  the  purpose  of  establishing  the 
reliability  of  various  analytical  methods  to  predict  performance 
and  to  provide  a  basis  for  improvements.  Such  studies  should 
emphasize  the  requirements  for  idealization  and  field  input 
information.  Without  this  information  analytical  techniques 
will  not  gain  acceptance  in  the  design  profession. 
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1. 
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Elasto-plas tic  analysis' of  a  thick-walled  circular  tube 
with  the  Von  Mises  yield  criterion. 

A  finite  element  idealization  of  the  thick-walled  circular 
tube  with  the  material  properties  used  in  the  analysis  is 
shown  in  Figure  1-1.  Typical  results  of  the  analysis  are 
presented  in  Figures  1-2  to  1-4. 


Elasto-plas  tic  analysis  of  a  circular  opening  with  the 
generalized  Mohr-Coulomb  yield  criterion. 
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with  the  material  properties  used  in  the  analysis  is 
shown  in  Figure  1-5.  Figure  1-6  shows  the  distribution 
of  vertical  and  horizontal  stresses  along  a  horizontal 
section  of  the  opening. 


3.  Combined  no  tension,  joint  perturnbation  and  elasto-plastic 
analysis  of  a  rectangular  underground  opening. 

A  finite  element  idealization  of  the  rectangular  opening 
with  the  material  properties  used  in  the  analysis  is  shown 
in  Figure  1-7.  Typical  results  of  the  analysis  are  presented 
in  Figures  1-8  and  1-9. 


APPENDIX  II  -  REFERENCES 

1.  Baker,  L.  E.,  Sandhu,  R.  S.,  and  Shieh,  W.  Y.,  "Application 
of  Elasto-plastic  Analysis  in  Rock  Mechanics  by  Finite 
Element  Method,"  Proc.  of  Eleventh  Symposium  on  Rock 
Mechanics,  Berkeley,  California,  June,  1969. 


hh 


2.  Brown,  Gilvert  L.,  Morgan,  Ernest  D. ,  and  Dodd,  Jerry  S. 

"Rock  Stabilization  at  Morrow  Point  Powerplant,  "Journal 
of  the  SM  and  FD,  ASCE,  Vol,  47,  No.  SMI,  Proc.  Paper 
7820,  January  1971. 

3.  Bureau  of  Reclamation,  "A  Technical  Representation  of  Mor¬ 
row  Point  Dam  and  Powerplant  Foundation  Investigation,  Colo¬ 
rado  River  Storage  Project,"  Denver,  Colorado,  1965. 

4.  Chang,  C.  Y.,  and  Nair,  K. ,  "A  Theoretical  Method  for  Evalu¬ 
ating  Stability  of  Openings  in  Rock,"  Final  Technical  Report 
for  Contract  Number  H0210046,  U.  S.  Bureau  of  Mines,  (May, 
1972). 

5  Clough,  G.  W.  and  Duncan,  J.  M.,  "Finite  Element  Analyses 
of  Port  Allen  and  Old  River  Locks,  "College  of  Engineering, 
Office  of  Research  Services,  University  of  California, 
Berkeley,  California  1969. 

6.  Dodd,  Jerry  S.,  "Morrow  Point  Underground  Powerplant  Rock 
Mechanics  Investigations,"  U.S.  Bureau  of  Reclamation 
Report,  March  1967. 

7.  Druckcr,  P.  C.  and  Prager,  W.,  "Soil  Mechanics  and  Plastic 
Analysis  or  Limit  Design,"  0.  Appl .  Math..  Vol.  10,  1952. 
pp.  157-165. 

8.  Goodman,  R.  ,  Taylor,  R.  ,  and  Brckke ,  T.,  "A  Model  for  the 
Mechanics  of  Jointed  Rock,"  Proceedings ,  American  Society 

of  Civil  Engineers,  Vol.  94,  No.  SM3,  May,  1968,  pp.  637-660. 

9.  Heuer,  R.  E.  and  Hendron,  A.  J.,  Jr.,  "Geomechanical  Model 
Study  of  the  Behavior  of  Underground  Openings  in  Rock  Sub¬ 
jected  to  Static  Loads,  Report  1,  Development  of  Modeling 
Techniques,  Oct.  1909,  Report  2,  Tests  on  Unlined  Openings 
in  Intact  Rock,  Feb.  1971,"  Dept,  of  Civ.  Engineering,  Uni¬ 
versity  of  Illinois,  Urbana,  Illinois. 

10.  Judd,  W.  R.  and  Perloff,  W.  H. ,  "Strain  Distribution  Around 
Underground  Openings  -  Tech.  Report  No.  5,  Comparison  Be¬ 
tween  Predicted  and  Measured  Displacements,"  Purdue  Univer¬ 
sity,  Sch.  of  Civ.  Engrg . ,  Lafayette,  Indiana,  March  1971. 

11.  Patton,  F.  D.  ,  "Multiple  Modes  of  Shear  Failure  in  Rock," 
Proc  1st  Congress  of  the  International  Society  of  Rock 
Mechanics,  Lisbon,  Vol.  I,  pp.  509-515,  1966. 

12.  Prager,  W.,  and  Hodge,  P.  G.,  "Theory  of  Perfectly  Plastic 
Solids,"  Published  by  John  Wiley  f,  Sons,  Inc.  1951. 


! 


*»5 


13.  Rcy  cs,  S.  F.  ,  "Bias t i c - Plas tic  Analysis  of  Underground 
Openings  by  the  Finite  Element  Method,"  thesis  presented 
to  the  University  of  Illinois,  Urbana,  in  1966,  in  partial 
fulfillment  of  the  requirements  for  the  degree  of  Doctor 
of  Phil sophy. 

14.  Reyes,  S.  F. ,  and  Deere,  D.  U.,  "Elastic-Plastic  Analysis 
of  Underground  Openings  by  the  Finite  Element  Method," 
Proceedings ,  First  Congress  of  the  International  Society 
of  Rock  Mechanics,  Lisbon,  Portugal,  1966. 

15.  Zienkicwicz,  0.  C.,  Valliappan,  S.,  and  King,  I.  P., 
"Stress  Analysis  of  Rock  as  a  'No  Tension'  Material," 
Geotechnique ,  Vol .  18,  March  1968. 

16.  Ziekicwicz,  0.  C.  ,  Valliappan,  S.,  and  King,  I.  P.  , 

"Elas to-Plastic  Solutions  of  Engineering  Problems: 

'Initial  Stress'  Finite  Element  Approach,"  Inter. 

Journal  for  Numerical  Methods  in  Engineering  (1969) 

Vol.  1,  No.  1,  75-100. 

17.  Ziekrcwicz,  0.  C.  ,  Best,  B.,  Dullage,  C.  and  Stagg,  K.  G. 
"Analysis  of  Nonlinear  Problems  in  Rock  Mechanics  with 
Particular  Reference  to  Jointed  Rock  Systems,"  Procs. 

2nd  Congress  of  the  International  Society  of  Rock 
Mechanics,  Vol.  3,  No.  8-14,  Belgrade,  Yugoslavia,  1970. 


O  KN  o 

o  •  o 

«  O  r-t  o 

II  II  II  II 

UJ  U  O 


■ 


lN3W3DV*ldS  1 0 
nviNoziaoH  on 


J 


FINITE  ELEMENT  MESH  FOR  AN  ELASTO-PLASTIC  ANALYSIS  OF  A 
THICK-WALLED  CIRCULAR  TUBE  (b  =  2a) 
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O  -  INCREMENT  2,  ITERATION  3,  p  =  1.00k 
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FIG  1-2  -  DISTRIBUTION  OF  CIRCUMFERENTIAL  STRESS 
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FIG.  1-5  -  FINITE  ELEMENT  MESH  FOR  AN  ELASTO- PLASTIC  ANALYSIS  OF 
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FIG  1-6  ■-  VERTICAL  AND  HORIZONTAL  STRESSES  ALONG  HORIZONTAL 

SECTION  OF  A  CIRCULAR  OPENING 


DEPTH  BELOW  SURFACE  (FT.)  NO  HORIZONTAL 

^-displacement 


-  — ? 


NCREMENT  PERCENTAGE  OF  TOTAL 


BOUNDARY  PRESSURE 


■NO  VERTICAL  DISPLACEMENT 
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APPENDIX  B 

Abstract  of  oral  presentation  "Roof  Movement  Due  to  Increasing 
Roof  Span  in  an  Underground  Opening"  given  at  the  Fourteenth  Symposium 
on  Rock  Mechanics,  held  at  the  Pennsylvania  State  University,  University 
Park,  PA,  June  11-14,  1972. 

This  presentation  included  results  of  research  performed  under 
contract  H0210017  "Stresses,  Deformations,  and  Progressive  Failure  in 
Non-Homogeneous  Rock"  with  Ohio  State  University,  Columbus,  Ohio,  Dr.  R.  S 
Sandhu,  Principal  Investigator.  Following  is  information  on  obtaining 
technical  reports  and  a  magnetic  tape  containing  the  computer  programs 
developed  under  this  contract  from  the  National  Technical  Information 

Service,  U.S.  Department  of  Commerce,  5265  Port  Royal  Road,  Springfield, 
Va.  22151. 
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Price 

Semiannual  Technical  Report 

AD-731  688 
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Final  Technical  Report  -  Vol.  1 

AD-768  648/8 

5.501 

Final  Technical  Report  -  Vol.  2 

AD-768  649/6 

4.501 

Final  Technical  Report  -  Vol.  3 
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2.751 

Magnetic  tape  containing  computer  programs 
dc-  eloped  under  this  contract 

AD-768  651/2 

250.00 
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$312.50 

Price  shown  is  for  paper  copy.  These  reports 
microfiche  at  $1.45  per  copy. 
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ROOF  MOVEMENT  DUE  TO  INCREASING  ROOF  SPAN  IN  AN  UNDERGROUND  OPENING 

by 

Dr.  l\.  Knrwoski1,  Dr.  R.  Sandhu2,  and  Dr.  R.  Singh3 

A  finite  element  analysis  is  made  of  the  roof  displacements  in  an  under¬ 
ground  room  as  the  roof  span  is  incrementally  increased. 

The  prototype  opening  was  sited  in  a  large  tertiary  gravel  deposit.  The 
material  was  described  as  cemented  gravel  with  a  bulk  density  of  130  lbs/ ft3. 
No  other  engineering  properties  were  available. 

Initially,  two  rooms  were  driven  off  an  adit.  The  rooms  were  driven 
side-by-side  on  40- foot  centers  to  a  length  of  75  feet.  One  room  was  pro- 
gresively  widened  from  10  feet  to  15,  22,  30,  and  40  feet  in  a  sequence  which 
increased  the  room  span  on  each  side  of  the  initial  room  centerline.  Each 
widening  was  done  along  the  entire  room  length  at  or  e  by  blasting.  This 
procedure  left  an  irregular  5-  to  10- foot  pillar  uetween  rooms. 

Roof  to  floor  closure  measurements  were  made  at  eight  stations  along 
the  initial  centerline  of  the  room.  In  total,  24  closure  measurements  were 
made. 

The  closure  measurements  show  the  classical  creep  curve  profile  between 
widening  events.  From  these  measurements,  the  "instantaneous"  sag  in  the  roof 
due  to  each  widen’ ng  event  was  deduced.  No  attempt  as  made  to  calculate 
the  creep  deformation  of  the  roof. 

Finite  elements  were  used  to  model  a  typical  prototype  cross-section  as 
a  2-D  plain  strain  condition  and  stressed  by  gravity  induced  loads.  The 
material  engineering  properties  used  in  the  model  were  estimated.  The  room 
widening  events  were  simulated  in  the  model  by  decreasing  the  stiffness  of 


2  Mining  Engineer.  Spokane  Mining  Research  Center. 

3  Associate  Professor,  Civil  Engineering,  Ohio  State  University. 
Formerly  Research  Associate,  Ohio  State  University. 
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elements  which  were  to  be  removed  by  the  excavation  sequence.  Results  were 
output  at  locations  in  the  model  to  compare  with  the  prototype  closure 
measurements. 

The  calculated  closures  compare  favorably  with  the  measured  values. 

It  is  felt  that  to  obtain  closer  comparison  improvements  must  be  made  in 
the  prototype  measurement  technique  and  in  the  material  property  estimates. 
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APPENDIX  C 

Manuscript  of  "A  Parametric  Study  of  Stresses  in  Steel  Support  for 
a  Tunnel"  presented  at  the  Fifteenth  Symposium  on  Rock  Mechanics  held 
at  the  South  Dakota  School  of  Mines ,  Rapid  City ,  South  DaKota , 

September  17-19,  1973.  This  presentation  included  results  of  research 
performed  under  contract  H0210017  "Stresses,  Deformations,  and  Progres¬ 
sive  Failure  in  Non-Homogeneous  Rock"  with  Ohio  State  University,  Columbus 
Ohio,  Dr.  R.  S.  Sandhu,  Principal  Investigator.  Information  on  obtaining 
reports  and  computer  programs  developed  under  this  contract  are  given 
on  page  55. 
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A  PARAMETRIC  STUDY  OF  STRESSES  IN  STEEL 
SUPPORT  FOR  A  TUNNEL 

by 

1  2 

i  Or.  Ranbir  S.  Sandhu  and  Dr.  William  J.  Karwoskl 
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ABSTRACT 

A  two-dimensional  finite  element  computer  code  is  applied  to  the  analysis  of 
stress  and  bending  moment  in  a  steel  arch  and  invert  strut  commonly  used  for 
(  ground  support.  Parametric  investigations  are  made  to  delineate  and  rank  rock 

physical  properties  that  affect  the  loads  imposed  on  the  ground  support  structure. 
In  the  analysis  it  is  assumed  that  elastic  rock  deformations  have  taken  place 
prior  to  support  installation.  The  support  is  loaded  by  assuming  that  zones  of 
"failed  rock"  develop  around  the  support  and  other  mining  loads  are  imposed  as 
boundary  loads  on  the  computational  mesh. 

The  effect  of  asymetrical  blocking  and  blocking  density  on  stresses  and  bend¬ 
ing  moments  in  the  steel  support  is  determined.  The  elastic  modulus  of  the  rock 
had  the  most  significant  effect  on  the  state  of  stress  in  the  support.  Increased 
blocking  density  reduces  bending  moment  stresses,  and  asymetrical  blocking 
Increases  bending  moment  stresses. 

'  Professor,  Department  of  Civil  Engineering,  Ohio  State  University,  Columbus, 
Ohio. 

\  ^  MTning  Engineer,  Spokane  Mining  Research  Center,  U.S.  Bureau  of  Mines,  Spokane, 

Wash. 
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INTRODUCTION 

The  purpose  of  this  investigation  is  to  determine  the  influence  of  rock  pro¬ 
perties  and  the  number  and  distribution  of  blocking  points  on  stresses  and  bending 

i 

moments  in  a  steel  arch  used  for  mine  support.  Young's  modulus,  Poisson's  ratio, 
cohesioft,  and  the  angle  of  internal  friction  are  varied  over  a  wide  range  of  vaiues 
to  study  the  stress  distribution  in  a  continuously  loaded  steel  arch.  Using  the 
combination  of  physical  properties  that  produce  the  maximum  bending  moment  and 
stresses,  three  cases  with  varying  point  load  distributions  on  steel  arches  are 
Investigated. 

Figure  1  shows  the  configuration  of  the  tunnel  opening  and  the  four  different 
blocking  systems  studied.  The  initial  stress  field  was  specified  as  hydrostatic 
pressure  corresponding  to  an  r>erburden  depth  of  1,000  feet  and  materiai  density 
of  165  pounds  per  cubic  foot.  The  range  of  materiai  property  values  is  shown  in 
table  i . 

All  data  on  the  initial  stress  field,  material  properties,  steel  arches, 
blocking  material,  and  locations  were  provided  by  the  U.S.  Bureau  of  Mines 
Spokane  Mining  Research  Center.  The  finite  element  computer  program  used  in  these 
Investigations  was  developed  under  ARPA  contract  number  H02100I7  with  Ohio  State 
University  Research  Foundation.  The  contract  final  report  (reference  1)  and  a  copy 
of  the  computer  program  may  be  obt??ned  from  the  Defense  Documentation  Center, 
Alexandria,  Virginia  223i*».  (See  p-ige  55) 

Rational  For  Model  Load  Development 

When  a  tunnel  Is  excavated,  the  load  carried  by  the  materiai  removed 
must  be  carried  by  the  rock  in  the  tunnel  walls  and  by  the  unexcavated  rock 
ahead  of  the  face.  Continued  excavation  at  the  face  results  in  "loss  of 
support"  further  increasing  the  stress  In  the  walls.  For  linear  homogeneous 
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Isotropic  elastic  rock,  It  has  been  shown  (reference  2)  that  theoretically  this 
effect  Is  felt  only  In  a  region  one  diameter  away  from  the  face.  If  supports 
are  Installed  Immediately  after  excavation,  they  will  share  this  transfer  of 

I 

load  as  the  face  is  advanced.  It  Is  observed  in  field  tests  (reference  2)  that 
rock  movement  continues  for  a  long  time  before  reaching  stabilization.  This 
results  In  continued  growth  in  the  load  transferred  to  the  supports.  Also, 
exposure  to  atmosphere,  loss  of  gouge  material,  and  blasting  damage  may  alter 
the  mechanical  properties  of  the  rock  mass,  resulting  In  Increasing  deformation 
and  Increasing  support  stresses.  The  load  development  may  be  associated 
with  one  or  more  of  the  following  mechanism*: 

a.  Upon  continued  excavation  at  the  face,  the  removal  of  rock  res  'ts 
In  Increased  rock  load  being  transferred  to  the  walls  of  the  portion  already 
excavated  and  the  supports  in  that  portion. 

b.  Time-dependent  deformation  of  rock  Is  resisted  by  the  supports 
resulting  In  their  taking  on  Increasing  load. 

c.  Change  In  material  properties  after  Installation  of  supports  will 
result  in  additional  deformation  which  in  turn  will  lead  to  Increased  stresses 
In  the  supports  resisting  such  deformation. 

d.  By  blocking,  a  prestress  may  be  introduced  to  support  rock.  Wedging 
of  the  blocks  will  give  equal  and  opposite  forces  acting  on  the  tunnel  surface 
and  the  support  structure. 

in  the  work  reported  herein,  the  influence  of  various  rock  properties 
(see  Table  I)  upon  support  stresses  was  studie4  assuming  load  development 
primarily  through  mechanicism  (a)  described.  The  system  In  figure  1 
(case  a)  was  analyzed  to  study  the  effect  of  variation  of  parameters  and 
to  rank  them  In  order  of  importance.  Calculations  were  made  varying  each 
of  the  parameters  between  the  limits  given  one  at  a  time.  The  "constant" 
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values  of  other  parameters  were  taken  as  the  midpoint  of  the  specified 
range.  The  combination  of  parameter  values  corresponding  to  the  worst 
stress  In  the  supports  was  then  used  to  analyze  cases  b,  c,  d  shown  In 
figure  1. 

ASSUMPTIONS  MADE  IN  THE  ANALYSIS 
Two  Dimensional  Finite  Model 

The  assumption  of  two-dimensional  symmetry  can  only  be  strictly  true  when 
the  tunnel  Is  continuously  supported,  for  example,  with  concrete  because  every 
plane  normal  section  through  the  tunnel  is  Identical.  When  the  tunnel  Is  supported 
Intermittently  with  steel  supports,  every  plane  normal  section  through  the 
tunnel  Is  not  Identical,  and  the  analysis  must  be  three-dimensional  to  be  truly 
representative.  Furthermore,  rock  reaction  to  line  loads  produced  by  steel  sets 
Is  dlstln'-  ly  different  to  uniform  loads  produced  by  a  continuous  lining.  Dixon 
(reference  3)  discusses  this  In  detail  and  reports  a  method  whereby  problems  of 
Intermittent  supports  In  tunnels  can  be  made  amenable  to  two-dimensional  analysis 
by  scaling  Young's  modulus  of  elasticity  for  the  rock  mass. 

In  this  paper  the  Interest  Is  In  the  trend  of  the  support  stresses  as  various 
rock  mass  properties  are  increased  or  decreased  In  value.  The  same  trends  are 
expected  to  exist  when  the  three-dimensional  aspects  of  the  problem  are  accurately 
modelled.  A  two-dimensional  model  of  the  tunnel  and  Its  support  Is  deemed  suitable 
for  these  preliminary  Investigations. 

Extent  of  the  Finite  Model 

When  an  underground  opening  Is  excavated,  changes  In  the  stress  field 
and  associated  deformations  occur  in  the  entire  rock  mass.  The  principle  of 
local  action  implies  that  these  changes  diminish  with  increasing  distance  from 
the  opening.  In  the  finite  element  model,  a  finite  region  Is  generally  consider¬ 
ed.  On  the  boundary  of  this  finite  model,  force  or  displacement  boundary 
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conditions  must  be  prescribed.  These  can  be  based  on  the  assumption  either 
of  no  change  In  the  stress  fieid  or  of  no  deformation.  Neither  of  the  two  is 
true  for  finite  distances  from  the  opening,  and  the  two  assumptions  in  fact  give 
bounds  to  the  correct  solution.  Nair  {reference  k)  and  Kuthawy  (reference  5) 
among  others,  have  studied  the  effect  of  lateral  dimension  of  the  finite  model 
and  of  the  choice  of  boundary  conditions  on  stresses  and  deformations 
In  the  vicinity  of  underground  openings.  For  the  present  study, 
allowing  for  the  specified  hydrostatic  initial  stress  fieid,  a  preliminary 
analysis  showed  that  it  would  be  adequate  to  model  a  region  extending  approxi¬ 
mately  seven  diameters  above  the  roof  of  the  opening,  five  diameters  horizon¬ 
tally  on  each  side  of  the  center-line  and  five  diameters  below  the  invert  of 
the  opening.  The  boundary  conditions  for  the  finite  element  model  are  Illus¬ 
trated  In  figure  2.  The  overburden  of  the  top  858  feet  was  replaced  by  an 
equivalent  vertical  load,  and  gravity  was  applied  to  the  mesh  Interior.  The 
stress  field  on  the  vertical  faces  AP,  CD  was  assumed  to  be  unaffected  by  the 
sequence  of  operations  In  the  opening  and  the  vertical  displacement  of  the  hori¬ 
zontal  section  BC  was  set  equal  to  zero. 

Modelling  of  Support  Structure  and  Sequence  of  Operations 

Tho  steel  supports  were  assumed  to  be  in  plane  stress  whereas  the  rock  and 
shotcrete  were  in  plane  strain.  The  spacing  of  steel  supports  was  specified  as 
three  feet.  Thus,  a  three-foot  length  of  the  tunnel  was  supported  by  each  support 
ring.  Figure  3  shows  a  typical  cross-section  used  in  the  finite  element 
model.  The  cross-section  of  the  steel  rib  was  represented  by  five  finite 
elements  to  obtain  reasonably  good  distribution  of  stresses  over  the 
cross-section.  The  shotcrete  and  the  steel  rib  were  assumed  to  be  bonded. 

If  there  is  no  bond  between  shotcrete  and  the  rib,  there  is  no  load  trans¬ 
ferred  through  sheai;  and  the  load  transfer  is  entirely  radial.  This  situa- 
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tion  would  be  similar  to  case  d,  in  figure  1,  except  that  the  blocking 

would  be  continuous.  For  cases  b,  c,  d,  the  wooden  blocks  were  assumed  to 

be  axial  members  not  transmitting  any  bending  moments. 

I 

The  shotcrete,  the  steel  rib,  and  the  timber  blocks  were  assigned  the 
properties  given  in  table  1.  It  was  assumed  that  C  and  $,  specified  for  the 
rock,  were  obtained  In  triaxiri  tests  on  cylindrical  specimens.  It  was 
assumed  that  the  rock  properties  do  not  change  as  a  result  of  excavation  of 
the  opening. 

1  Starting  with  the  initial  stress  state,  the  sequence  of  opening  excavation, 
installation  of  support,  and  support  ioad  development  is  simulated  in  the  compu¬ 
tation.  As  opening  excavation  is  made,  all  deformation  due  to  the  interaction 
between  the  opening  and  the  initial  stress  state  takes  place  at  that  time.  The 
steel  support  Is  placed  in  the  deformed  opening.  At  this  point  in  the  computa¬ 
tion  no  loads  would  develop  on  the  support  unless  the  existing  stress  state  Is 
changed.  This  is  done  by  Increasing  the  density  of  the  rock  surrounding  the 
support  from  85.2  ibs/ft^,  used  in  the  initial  computation,  to  i 65  ibs/ft^,  used 
in  the  computation  with  supports.  The  shaded  areas  in  figure  k  are  the  regions 
In  which  this  weight  density  change  was  made.  This  sequence  of  operations  is 
rationalized  as  simulation  of  model  ioad  development  mechanism  (a)  described 
earlier.  The  idea  that  the  weight  supported  by  the  steel  sets  is  limited  to  a 

l 

finite  region  of  rock  around  the  tunnel  is  taken  from  Kar'  Terzaghl's  "Rock 
Defects  and  Loads  on  Tunnel  Supports"  (Proctor  and  White,  reference  6).  The 
geometry  of  the  shaded  areas  in  figure  4  corresponds  to  some  observed  tunnel  pro¬ 
files  where  the  rock  has  fallen  out  of  the  roof.  Because  good  information  on 
how  rock  loads  develop  is  not  available,  it  was  arbitrarily  assumed  that  the  dif¬ 
ference  in  excavation  profiles  in  figure  i  case  b  and  case  a  represented  the 
rock  load  which  could  develop  for  cases  a  and  d.  For  case  b,  the  rock  exertlnr 
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the  gravity  load  Is  shown  shaded  In  figure  *»b,  the  extent  was  arbitrarily  taken 
as  an  average  of  3*5  feet  thickness  beyond  the  excavation  line.  Similarly  for 
case  c,  the  extent  was  arbitrarily  taken  as  an  average  of  3.5  feet  thickness 
beyond  the  excavation  line. 

A  number  of  Important,  and  perhaps  novel  assumptions,  was  made  In  order 
to  satisfy  the  problem  Input  requirements.  This  Is  a  particularly  difficult 
task  especially  In  areas  such  as  rock  load  transfer  or  rock  load  development 
where  knowledge  Is  scarce.  However,  In  view  of  the  moderate  goals  of  the  research 
work,  this  Is  acceptable. 

RESULTS 

Prel Imlnary 

In  studying  the  Importance  of  material  parameters,  the  quantities  of 
Interest  were  the  rru.  Imum  stresses  In  the  steel  member.  It  Is  customary 
to  study  the  axial  and  shearing  forces  and  bending  moments  In  such  structural 
components.  Hence,  these  quantities  were  determined.  It  should  be  noted 
that  the  basic  output  from  the  finite  element  analysis  Is  the  components  of 
stress  evaluated  at  the  center  of  each  of  the  five  elements  Into  which  the 
steel  member  is  divided.  Tht  moments  and  forces  were  obtained  by  numerical 
tntearatton  of  the  stress  values. 

Ftgure  5  shows  the  cross-sections  at  which  the  stresses  as  well  as  the 
forces  and  moments  In  the  steel  rib  were  computed. 

Influence  of  Material  Properties 

Four  material  parameters,  viz.,  Young's  modulus,  Poisson's  ratio,  cohesion 
and  angle  of  Internal  friction  were  to  be  considered  to  establish  their  order  of 
Importance  In  terms  of  their  Influence  on  the  support  forces. 
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For  the  specified  range  of  the  rock  properties,  the  cohesion  and  the  angle 
of  Internal  friction  did  not  Influence  the  stresses  In  the  structural  supports. 
Figure  6  shows  a  plot  of  Jj  ■  Oj  +  +  0^,  ^rst  Invariant  of  the  st"ess 

tensor  as  the  abclssa  and  — *.  /g2.  ?)) — t-J?) - £lL  ]  ^  ,  the  second 

Invariant  of  the  stress  deviator  tensor  as  the  ordinate.  Generalized  Mohr-Cou- 
lomb  yield  criteria  corresponding  to  the  prescribed  range  of  values  of  cohesion 
and  angle  of  Internal  friction  are  shown  and  also  the  stress  paths  traced  by 
points  of  critical  locations  around  the  tunnel  opening. 

The  generalized  Mohr-Coulomb  yield  law  Is, 


i 

a  Ji  +  J2  -  k 


,  1  2  s  I  n  4> 

where:  a  ■  rs  -  siV) 

,  1  6C  cos  4> 

k  "  n  13  -  sin#)' 


<$>  angle  of  Internal  friction 
C  ■  cohesion 

C  and  ^  are  determined  from  a  trlaxlal  test  on  a  cylindrical  specimen  under 
axlsymmetric  radial  stress.  Derivation  of  the  above  equations  Is  given  by  Singh 
(reference  5). 

In  Figure  6,  paths  A,  B,  C,  D,  E,  are  traced  by  elements  located  around  the 
face  of  the  underground  opening.  Points  A,  C  refer  to  elements  at  the  Invert  and 
the  crown,  respectively,  and  points  B,  D,  E  correspond  to  elements  on  the  side  of 
the  opening.  The  locations  are  Indicated  In  figure  5.  The  Initial  state  in  ail 
cases  Is  of  hydrostatic  stress.  The  terminal  points  represent  the  state  after 
excavation.  The  development  of  rock  load  has  little  influence  on  the  stress  state 
In  rock  for  the  specified  values  of  Young's  modulus  and  Poisson's  ratio.  For  all 
locations,  the  entire  stress  history  was  found  to  be  well  below  the  yield  criteria 
corresponding  to  the  range  of  parameter  values  specified  in  Table  I. 
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Figures  7(a)  through  7(c)  show  the  Influence  of  variation  of  Young's  modulus 
upon  bending  moments  and  upon  axial  and  shearing  forces.  Only  sections  with  the 
worst  forces  have  been  plotted.  The  maximum  bending  moments  were  at  section  16 
(crown),  the  maximum  axial  force  at  section  5  (side),  and  the  maximum  shearing  force 
at  section  10.  In  all  cases,  decrease  In  elastic  modulus  of  rock  was  seen  to 
result  In  Increased  moments  and  forces.  This  Is  Indeed  to  be  expected.  The 
decrease  In  support  forces,  with  Increasing  Young's  modulus,  Is  rapid  at  first 
and  then  Is  less  pronounced.  This  Is  due  to  the  fact  that  rock  «*.'  -‘ormatlon  Is 
proportional  to  reciprocal  of  the  modulus.  Also,  for  very  large  moduli,  the 
strains  are  extremely  small, and  difficulties  arise  with  computational  precision. 

Figures  8(a)  through  8(c)  show  the  influence  of  the  variation  in  the  Poisson's 
ratio  on  the  moments  and  axial  and  shear  forces  in  the  supports.  Again,  only 
results  for  the  worst  sections  have  been  plotted.  In  the  crown  and  the  side 
(section  16  and  5),  an  Increase  in  Poisson's  ratio  results  in  decreased  bend¬ 
ing  moments  whereas  at  section  10  the  bending  moment  Increases  somewhat. 

For  tne  specified  range  of  values  for  the  material  parameters,  the  bending 
moments  and  the  shearing  forces  were  small.  The  major  effect  was  the 
•xtal  force  In  the  member.  An  explanation  for  the  bending  effect  being  very 
small  may  be  found  in  the  assumption  that  the  steel  member  is  restrained  by 
the  shotcrete.  Significant  bending  of  the  steel  member  must  involve  significant 
Changes  In  curvature.  This  is  not  possible  for  a  member  restrained  from  radial 
movement  by  relatively  unyielding  rock.  Also,  the  rock  properties  are  assumed 
to  be  unaffected  by  the  excavation  process.  This  may  not  be  true.  There  can  be 
considerable  change  in  deformabi  1  i ty  of  a  rock  mass  due  to  excavation  of  the 
underground  opening. 
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In  this  study,  It  ts  evident  that  supports  In  lower  elastic  modulus  rock 
will  develop  greater  stresses.  High  value  of  Poisson's  ratio  has  the  effect  of 
redistributing  stresses  decreasing  the  peaks  and  Increasing  the  lowest  values. 

Using  the  lowest  value,  in  the  prescribed  range,  of  Young's  modulus  and  the 
mlrtmum  as  well  as  the  maximum  values  of  Poisson's  ratio,  stresses  for  cases  a,  b, 
c,  d  of  different  excavation  profiles  and  blocking  arrangements  were  determined. 
The  resulting  stress  distributions  are  shown  in  figure  9. 

Additional  Studies 

Additional  studies  using  values  of  Young's  modulus  less  than  the  minimum 
of  the  range  given  In  table  1  were  carried  out  to  consider  situations  where  the 
elasttc  modulus  of  rock  In  the  vicinity  of  an  underground  opening  may  be  signifi¬ 
cantly  reduced  by  damage  during  excavation  or  deterioration  with  exposure  over  a 
long  time.  The  simulation  assumed  that  after  installation  of  supports  and  place¬ 
ment  of  shotcrete,  Young's  modulus  may  reduce  from  an  Initial  value  of  1  x  10^ 
pounds  per  square  Inch.  Different  terminal  values  of  Young's  modulus  used  were 
0.75  x  10  ,  0.67  x  10  ,  0.4  x  10  ,  0.25  x  10°  pounds  per  square  inch.  Poisson's 
ratio  was  assumed  .o  be  0.49  throughout.  For  case  a,  the  results  are  plotted 
In  figures  10(a;  through  10(c).  As  might  be  expected,  greater  reduction  In  the 
elasttc  modulus  was  associated  with  greater  bending  moments  and  axial  and  shear 
forces.  For  reduction  of  Young's  modulus  to  0.25  x  10^  psl,  the  maximum  longitu¬ 
dinal  stress  would  be  over  20,000  psl.  Considering  that  deterioration  of  rock  Is 
more  likely  to  occur  in  cases  b  and  d,  analysis  was  performed  corresponding  to 
a  reduction  In  Young's  modulus  from  1  x  10^  to  0.4  x  10^  pounds  per  square  Inch. 
Figure  11  shows  the  distribution  of  longitudinal  stresses  at  critical  sections 
for  the  cases  a,  b  and  d. 
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CONCLUSIONS 

The  tendency  of  the  rock  to  deform  Is  a  key  factor  In  load  transfer  to  the 
structural  support  system.  Higher  modulus  Implies  less  deformation  for  the 
same  rock  load  and  consequently  less  load  transferred  to  the  steel  member. 

Higher  values,  tending  towards  .5,  of  Poisson's  ratio  are  associated  with 
redistribution  of  stress  In  the  rock.  This  stress  redistribution  Is  reflected 
In  a  more  uniform  stress  distribution  In  the  steel  support. 

Sparce  blocking  Increases  bending  moments  In  the  support  structure.  The 
total  load  transferred  to  the  support  structure  decreases  as  the  amount  of  block 
Ing  between  the  steel  set  and  rock  Increases.  The  wood  blocking  acts  much  like 
a  "back  packing"  material. 

The  results  of  these  analyses  appear  to  be  Intuitively  correct.  The  next 
step  Is  to  use  the  computer  code  to  analyze  a  documented  field  case. 
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FIGURE  1.  -  Tunnel  Excavation  and  Support  Blocking  Configuration 
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FIGURE  2.  -  Boundary  Conditions  on  the  Finite  Element  Model. 
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(a)  longitudinal  Cross-Section  of  Tunnel. 
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(b)  Representative  Slice. 


FIGURE  3.  ”  Longitudinal  Cross-Section  of  runnel  and  Representative 

Slice. 


(a)  Tunnel  Cross-Sect  ton  Showing  Rock  Load  for  Cases  a  and  d. 


(b)  Tunnel  Cross-Section  Showing  Rock  Load  for  Case  b. 


(c)  Tunnel  Cross-Section  Showing  Rock  Load  for  Case  c. 


FIGURE  A.  -  Rock  Load  on  Supports. 


FIGURE  5.  ~  Location  of  Structure  Sections 
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FIGURE  6.  -  Stress  Paths  for  Points  Around  the  Underground  Opening 
(for  location  of  points  see  FIG.  5). 
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FIGURE  $.  -Stress  Distribution  For  Critical  Combination  of  Rock  Properties. 


Section  25  Section  25 


PARAMETERS 

MATERIAL 

(a)  Rock 

(b)  Shotcrete 

(c)  Steel 

Young's  Modulus, 

E  (psl) 

IxlO6  -  10xl06 

2.71  x  106 

30  x  106 

1.5  x  106 

’olsson's  Ratio, 

V 

0.1  -  0.5 

0.1 

0.231 

0.03 

Cohesion,  c  (psl) 

1000  -  5000 

3300 

— 

— 

Angle  of  Internal 

Fr let  Ion,  (tan  $ ) 

1.0  -  2.1 

1.0 

— 

— 

later  la  1  Density 
(pcf) 

165.0 

150.0 

A90.0 

27.0 

(a)  Data  furnished  by  sponsor. 


(b)  "Manual  of  Concrete  Practice,"  Part  2,  1968,  American  Concrete 
institute 

E  -  w1,5  33  /f ' 

where  w  -  150  Ib./c.ft. 

f'  ■  compressive  strength,  assumed  as  2000  psl 
c 

(c)  "Manual  of  Steel  Construction,"  AISC 

(d)  "Wood  Handbook,"  Forest  Products  Laboratory,  Forest  Research. 
The  properties  are  for  western  white  pine. 


Table  1.  -  Material  Properties. 


83 


APPENDIX  D 

Summary  of  finite  element  computer  code  developed  under  ARPA  contract 
H0220035  "Analytic  Modeling  of  Rock-Structure  Interaction"  with  Agbabian 
Associates,  El  Segundo,  Calif.,  Dr.  Jeremy  Isenberg,  Principal  Investiga¬ 
tor.  The  computer  program  and  all  technical  reports  generated  under 
this  contract  may  be  obtained  from  the  National  Technical  Information 
Service,  U.S.  Department  of  Commerce,  5285  Port  Royal  Road,  Springfield, 
Va.  22151.  Inforaation  required  for  ordering  these  materials  are  given 
below. 


AD  No.  Price 


Semiannual  Technical  Report 

AD-749 

373 

$  3.00* 

Annual  Technical  Report  -  Vol.  1 
(describes  range  of  possible  application  and 
discusses  aspects  of  the  program  and  its 
theoretical  basis) 

AD-761 

648 

6.751 

Annual  Technical  Report  -  Vol.  2 
(explains  the  overall  operation  of  the 
program) 

AD-761 

649 

5.451 

Annual  Technical  Report  -  Vol.  3 

AD-761 

650 

3.001 

Magnetic  tape  containing  the  computer  program 

AD-762 

422 

97.00 

Foreign 

$122.50 

1 Price  shown  is  for  paper  copy,  microfiche  copies  are  also  available  at 
$1.45  per  report. 
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APPENDIX  D 

Summary  of  Finite  Element  Computer  Code 

Existing  finite  element  technology  was  consolidated  into  a  single,  gen¬ 
eral  purpose  computer  program.  This  computer  program  considers  static  two- 
and  three-dimensional,  nonlinear  analysis.  Gravity  loading  and  sequence  of 
excavation  of  construction  are  also  considered.  The  following  types  of  ele¬ 
ments  are  available,  (a)  beam,  (b)  rod,  (c)  axisymetric  quad,  (d)  plane 
and,  (e)  three-dimensional  hexadedron,  (f)  thick  shell,  (g)  plane  slip  or 
Joint,  (h)  axisymetric  slip  or  joint. 

All  elements  except  the  beam  and  thick  shell  have  nonlinear  properties 
which  may  be  assigned  are:  (a)  variable  moduli,  (b)  elastic  or  variable  modu¬ 
lus  with  ideal  plasticity.  Yield  function  may  be  isotropic  or  anisotropic, 
(c)  nonlinear  elastic  with  ideally  plastic  fracture  with  anisotropic  fracture 
criterion,  (e)  viscoelasticity  (Kelvin,  Maxwell,  three-parameter  fluid), 

(f)  viscoplasticity. 

The  computer  program  was  intended  to  solve  problems  in  rock  mechanics 
and  mining  engineering.  Thus,  the  excavation  and  construction  options  are 
intended  to  aid  in  analysis  of  ehafts  or  chambers  in  rock  having  nonlinear 
properties  from  the  start  of  construction  through  excavation  and  installa¬ 
tion  of  temporary  supports  to  the  installation  of  permanent  supports. 

The  program  uses  the  direct  stiffness  method  of  structural  analysis. 
Degrees  of  freedom  are  defined  at  nodal  points  which  are  at  corners  of  ele¬ 
ments.  Element  stiffnesses  are  deposited  in  a  global  stiffness  matrix. 

Nodal  point  forces,  which  include  external  loads  and  internal  resisting 
forces  are  expressed  in  a  global  load  vector.  A  form  of  Chloeski  decom¬ 
position,  involving  triangularization  of  the  global  stiffness  matrix,  reduc- 
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tion  of  the  load  vector  and  back  substitution  of  the  reduced  load  vector 
Into  the  trlangularized  stiffness  matrix,  Is  used  to  obtain  Incremental 
displacements.  The  program  has  several  options  for  performing  these 
operations  from  which  it  chooses  depending  on  the  available  core.  According 
to  one  option,  the  stiffness  matrix  is  formed  block  by  block  (one  bandwidth 
wide)  and  stored  on  peripheral  units;  it  is  then  retrieved  block  by  block 
and  processed.  Under  a  second  option,  the  stiffness  matrix  and  load  vector 
are  formed  and  reduced  in  a  block-by-block  manner  and  written  on  peripheral 
storage  in  a  reduced  form.  Part  of  the  reduction  is  performed  in  local 
double  precision. 

These  operations  require  many  accesses  to  peripheral  storage.  The  number 
of  accesses  and  hence  the  computer  time  required  for  large  problems  would 
become  prohibitive  on  small  or  medium  sized  computers.  In  the  present 
program,  a  multituffering  technique  is  used  which  allows  data  to  be  trans¬ 
ferred  to  and  from  peripheral  storage  while  computations  are  being  per¬ 
formed  In  the  main  core  area.  This  procedure  virtually  eliminates  wait 
time  for  Input/Output  operations. 

The  computer  program  was  applied  to  the  analysis  of  floor  and  roof 
movements  during  the  construction  of  an  underground  hoist  room.  Good 
comparisons  were  made  between  computed  and  field  measurements. 


